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Abstract 

We present the first inelastic neutron scattering study of tiie sliort wavelengtii dynamics in a pliospiiolipid bilayer. 
We siiow tliat inelastic neutron scattering using a triple-axis spectrometer at the high flux reactor of the ILL yields 
the necessary resolution and signal to determine the dynamics of model membranes. The results can quantitatively 
be compared to recent Molecular Dynamics simulations. Reflectivity, in-plane correlations and the corresponding 
dynamics can be measured simultaneously to gain a maximum amount of information. With this method, dispersion 
relations can be measured with a high energy resolution. Structure and dynamics in phospholipid bilayers, and 
the relation between them, can be studied on a molecular length scale. 
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Motions in biological membranes span a wide 
range of length and time scales [1]. Short wave- 
length density fluctuations in the plane of the 
bilayer (as opposed to bilayer bending modes) 
are likely to play a key role in the transport of 
small molecules across the bilayer [2], but have 
received little experimental attention. While the 
structure of lipid bilayers has been the object 
of several studies in the last decades, dynamical 
properties even of simple model systems such as 
l,2-dimyristoyl-sn-glycero-3-phoshatidylcholine 
(DMPC) remain largely unknown. A more quan- 
titative understanding of short range molecular 
motions in a lipid bilayer is of fundamental interest 



in membrane biophysics. Experimental results can 
be linked to modern Molecular Dynamics (MD) 
simulations [3]. In this context, Chen et al. have 
recently presented a seminal measurement of the 
dispersion relation of the acyl chain fluctuations 
in lipids [4] by inelastic X-ray scattering (IXS). As 
an important result, Chen et al. found a minimum 
in the dispersion relation at about Qo=1.4 A~^, 
at the position of the static structure factor max- 
imum S(Q2;y). This soft mode was hypothetically 
linked to transport phenomena across the bilayer. 
Inelastic neutron scattering using a triple-axis 
spectrometer should offer a much better resolu- 
tion and avoid the danger of radiation damages 



Preprint submitted to Elsevier Science 



2 February 2008 



Analyzer 




Fig. 1. Accessible {Q,f2) range for a fixed energy of the 
scattered beam of Ej=10 meV. In the insert, a schematic 
of a triple-axis spectrometer is shown. 

in the sample, but has to date not been apphed 
to study the dynamics of hpid bilayers. A further 
important advantage of INS over IXS is the fact 
that different parts of the bilayer can be probed 
by use of selective deuteration. We carried out the 
first inelastic neutron scattering experiments in 
model membranes employing the cold triple-axis 
spectrometer IN12 at the high flux reactor at the 
ILL. In this study we used fully deuterated chains 
DMPC -d54 (Avanti hpids, Alabama) in D2O to 
minimize incoherent background. We report on 
experimental details and show examples of neu- 
tron diffraction and inelastic neutron scattering in 
deuterated DMPC bilayers. 

From a solution of DMPC -d54 (deuterated 
1 , 2-dimyristoyl-sn-glycero-3-phoshatidylcholine) 
in TFE/chloroform (1:1) with a concentration 
of 40 mg/ml, we prepared highly oriented mem- 
brane stacks on 2" silicon wafers (thickness 300 
fim) [5] . Ten of these wafers were stacked on top of 
and aligned with respect to each other to create 
a 'sandwich sample' consisting of several thou- 
sands of highly oriented lipid bilayers with a total 
mass of about 300 mg of deuterated DMPC. Alu- 
minum spacers between the substrates allow to 
let the heavy water vapor in between the wafers 
to hydrate the bilayers. In rocking scans we find a 
mosaicity of the sandwich of Aa;=0.6 deg. 

The insert in Fig. 1 shows a schematic of a triple- 
axis spectrometer. By varying the three axes of the 
instrument, the axes of rotation of the monochro- 



mator, the sample and the analyzer, it is possible to 
change the wavevectors k; and k f and the energies 
Ei and E/ of the incident and the scattered beam, 
respectively. The accessible (Q,ri) range for a fixed 
energy of the scattered beam E/ of 10 meV is shown 
in Fig. 1. It is limited by the range of incident neu- 
tron energies offered by the neutron guide as well 
as by mechanical restrictions of the spectrometer. 
The instrumental energy resolution in this config- 
uration is Ail=500 /icV. By choosing smaller inci- 
dent energies and energy transfers the energy reso- 
lution can be enhanced. By rotating the sample the 
scattering vector Q can be placed either within the 
plane of the membranes (Qxy) to measure inter- 
chain correlations or perpendicular to it (Q2) to 
measure reflectivity (see insert in Fig. 4). The use 
of a triple-axis spectrometer offers the unique pos- 
sibility to measure reflectivity, the static structure 
factor S(Qxy) in the plane of the membranes and 
in plane dynamics {S{Qxy, i^)) on the same instru- 
ment in the same run. To control temperature and 
the degree of hydration during the experiment the 
sample was kept in a humidity chamber. Inside the 
chamber there was a heavy water reservoir, ther- 
mally isolated against the sample. By adjusting the 
temperature of the chamber and independently the 
temperature of the reservoir using two Haake bath 
controllers, the humidity inside can be controlled. 
Close to the sample there were a humidity sensor 
and a platinum resistor installed. IN12 is equipped 
with a vertically focusing graphite monochromator 
and a vertically and horizontally curved graphite 
analyzer. The flux at the sample position is of the 
order of 5T0^ n/(cm^- s). For the inelastic mea- 
surements the counting times were about five min- 
utes per point. Figure 2 shows an energy-scan at 
Qxy=1.0 and T=30 °C, in the fluid phase of 
the DMPC bilayer. The signal consists of a sharp 
central peak and a broad contribution, centered at 
rt—0 meV. The solid lines are guides to the eye. At 
ri=-5.5 meV, the inelastic signal of the membrane 
is observed in form of a pronounced peak with a 
width (FWHM) of about Af7=l mcV. Figure 3 
shows the corresponding dispersion relation in the 
fluid La phase at T=30 °C as measured by several 
constant-Q scans. The measurements can quanti- 
tatively be compared to corresponding Molecular 
Dynamics simulations by Tarek et al. [3] . Figure 4 
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Fig. 2. Energy scan at Q=1.0 and T=30 °C. The insert 
shows the excitation of the DMPC-layers in magnification. 
Solid lines are guides to the eye. 
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Fig. 3. Dispersion relation in the fluid phase of the lipid 
bilayer. Note the minimum at the maximum of the static 
structure factor S(Q3;j/) shown in Fig. 4. 
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Fig. 4. The inter-acyl-chain correlation peak (left) and re- 
flectivity (right). The insert schematically shows the ori- 
entation of the sample with respect to the incoming beam 
for the in-plane measurements. Solid lines are fits after 
the Lorentzian- (chain-peak) and Gaussian-model (reflec- 
tivity). 

be modelled. At T=30 °C we find a next-neighbor 
distance of the acyl chains of da;j,=4.5 A and a cor- 
relation length of ^a;i/=ll A. From the reflectivity 
there is a distance between the bilayers of d2=55 A. 
The combination of diffraction and inelastic mea- 
surements leads to a complete picture of structure 
and dynamics of model membranes on a molecular 
length scale. This will be especially useful in tem- 
perature dependent measurements in the range of 
the gel-liquid phase transition to study the devel- 
opment of structure and dynamics and the relation 
between them. This work has been funded by the 
German Research Ministry under contract number 
05-300-CJB-6 (MCR). 



shows the corresponding inter-acyl-chain peak for 
temperatures below and above the gel-liquid phase 
transition (phase transition temperature Tc=21.5 
°C) and the reflectivity at T=30 °C. When heating 
from the more rigid gel-phase into the fluid phase, 
the inter-chain peak broadens drastically (decreas- 
ing correlation length S^^^y) and the peak center 
shifts to smaller Q^^y-values (larger average next 
neighbor distance and smaller packing density). 
The reflectivity shows five well-developed Bragg 
peaks giving information about the distance of the 
layers in z-direction. From the structure factors of 
the Bragg peaks, the z-profile of the bilayer can 
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